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Abstract 

Predictions for the out-going muon fluxes from the annihilation of neutralinos in 
, the center of the Sun and the Earth in supergravity unification are given. Effects of 

> ; 

uncertainties of the input data are analysed. It is shown that the out-going muon 

a\ \ 

\ flux measurements from the Sun and the Earth are complementary, with the Earth 

o 

o\ 



providing a larger flux for low values of fine tuning and the Sun providing a larger flux 
for high values of fine tuning, and that only a combination of the two can compete 
favorably with the direct detection measurements. The predictions are compared 



^ ; with the recent hmits from BAIKAL NT-96, MACRO, and BAKSAN. 



Recently analyses of the data from several neutrino telescopes (BAIKAL NT- 
960,MACRoi, BAKSAnI) have produced bench mark limits on the out-going muon fluxes 
from the annihilation of dark matter wimps in the centers of the Sun and the EarthH. 
In this Letter we give predictions for these out-going muon fluxes in supergravity grand 
unificationi'l taking into account the uncertainties in the local wimp density as well as the 
uncertainties in the rms wimp velocity and we compare these predictions with the recent 
experimental limitsEFl. We study the relative merit of using the Sun vs the Earth for the de- 
tection of dark matter as the SUSY spectrum gets progressively heavier and the flne tuning 
parameter of the radiative electro-weak symmetry breaking gets progressively larger. We 
also make a similar comparison of the direct vs the indirect detection of neutralino dark 
matter and show that only a combination of the out-going flux measurements from the Sun 
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and the Earth can compete effectively with the direct detection of neutrahno dark matter 
in mapping the parameter space of minimal supergravity. 

Supergravity unified models are extensions of the Standard Model, consistent with all 
current data, and under the constraint of R partly invariance lead to the existence of a cos- 
mologically stable neutralino x5 which can act as a candidate for cold dark matter. In these 
models the breaking of supersymmetry takes place in the hidden sector and the breaking 
is communicated via gravitational interactions to the visible sector where the breaking of 
supersymmetry can be parametrized by the following set for the minimal case (mSUGRA)i: 
the universal scalar mass (mo) at the GUT scale, the universal gaugino mass {mi/2) at the 
GUT scale, the universal trilinear coupling Aq at the GUT scale, tanP =< H2 > / < Hi >, 
where < H2 > gives mass to the up quarks and < Hi > gives mass to the down quarks and 
the leptons, and sign(/i) where fiHiH2 is the Higgs mixing term in the superpotential. De- 
tailed relic density analyses show that supergravity models can produce a wimp relic density 
consistent with a range of current cosmological models. We shall assume this relic density 
to lie in the range 0.05 < ^j^oh"^ < 0.3, where fl^o = p^o/pc, with p^o the average density 
of neutralinos in the universe and pc = SH"^ /SttGn the critical matter density, where Gn is 
the Newtonian constant and h is the Hubble parameter in units of 100kmsec~^ Mpc~^ . 

The basic elements of the analysis of the physics of out-going muons are the following: 
The Milky Way neutralinos are trapped by scattering and gravitational pull by the Earth 
and the Sun and accumulate in their cores where they eventually annihilate. Some of the 
final annihilation remnants contain neutrinos which propagate and undergo charge current 
interactions in the rock surrounding the detector and produce out-going muons which are 
detected. The basic theory of the various elements of the analysis are fairly well laid out 
and we follow the standard literature hereSil. However, we point out some basic features 
where specific elements of the SUGRA models enter. The flux of the out-going muons is 
given by $^ = F^/, where Ta is the rate of Xi ~ Xi annihilation in the center of the Earth 
or the Sun and f is the product of remaining factorsS^. A significant portion of the SUSY 
dependence in contained in Ta which is given analytically byil 
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= ^tanh\t/T) (1) 

Here C is the capture rate of the neutrahnos which depends on the wimp rehc density and 
on the wimp-nucleus cross section, t is the age of the Earth or the Sun and r = [CCa]'^^'^ 
where Ca depends on the wimp annihilation cross-section. In the limit t » t, one finds 

~ C/2, i.e., an equilibrium is reached between annihilation and capture. The condition 
of equilibrium depends on the particulars of the SUSY parameter space. For the case of the 
Sun one finds that the equilibrium condition is realized over essentially all of the parameter 
space of the model while for the case of the Earth this condition is highly model dependent. 

The predictions of the out-going muon fluxes in supergravity unifed models will be con- 
strained by several factors. In addition to the radiative electro-weak symmetry breakingi 
and the relic density constraints we impose the constraint from the flavor changing neu- 
tral current process 6 — > s + 7. This process is known to constrain dark matter analyses 
severely by eliminating large segments of the mSUGRA parameter space via the relic den- 
sity constraint^. The parameter space of mSUGRA will be further constrained by the 
naturalness criterion which we take to imply mo < 1 TeV, rrig < 1 TeV, tanjS < 25, and 
—7 < {At/mo) < 7, where At is the value of Aq at the electro-weak scale for the top 
quark channel. In the analysis of the relic density we use the accurate method which in- 
cludes the correct thermal averaging over the Breit-Wigner poles in the neutralino s channel 
annihilation^. 

In the analysis below we study the dependence of our predictions for the out-going 
muon fluxes from the Earth and the Sun including the uncertainties of the neutralino relic 
density in the Milky Way. Current estimates show that the local wimp density can vary in 
the rang (0.2 — 0.7)GeVcm ^. We parametrize this variation by ,^ = p^o/po which for 
Po = O.SGeVcm"^ lies in the range 0.7 < ^ < 2.3. Further, we study the dependence of the 
out-going muon fluxes on the uncertainties in the wimp velocitylli the details of which are 
given below. There could also be effects on the out-going muon fluxes due to rotation of the 
galaxy0. However, such effects may be small (e.g. an analysis of these effect on the direct 



detection rates of dark matter was found to be only ~ 10%). A brief discussion of the effects 
on the prediction of the muon fluxes under the constraint of the annual modulation signal 
claimed by the DAMA Collaborationlii is also given. Finally, we give an analysis of the 
out-going muon flux from the Sun and the Earth as a function of the fine tuning parameter 
and discuss the relative merits of the direct vs the indirect detection. In the analysis below 
we neglect the v oscillation effects. These effects could modify the detection rates from 
10% — 50% but are highly model dependent on the nature of the assumed v oscillation0. 
Inclusion of these effects would not change the conclusions of this work. 

We discuss now the results of our analysis. In Fig. (la) a plot of the maximum out-going 
muon flux for the case of the Earth is given as a function of m^o for three ^ values: ^ = 0.7, 
^ = 1 and ^ = 2.3. One finds that the maximum out-going muon flux is generally peaked for 
m^o below 60 GeV after which it shows a fall off with increasing vti^^. Interestingly one finds 
that the MACRO limits are already stringent enough to put constraints on mSUGRA in this 
case although the part of the parameter space so constrained is rather small. In Fig. (lb) a 
similar plot is given for the case of the out-going muon flux from the Sun. However, here the 
MACRO limits do not put any constraint on mSUGRA. In fact, the experimental sensitivity 
will have to increase by a factor of 10 or more before mSUGRA can be constrained by the 
Sun flux. A comparison of Fig. (la) and Fig. (lb) shows that for relatively low values of mj^o 
the Earth provides a better source for wimp detection relative to the Sun. We shall return 
to a more detailed comparison of the goodness of the Earth vs the Sun for wimp detection 
in the discussion of Fig. 3. 

The analysis of Figs. (la) and Fig. (lb) was based on a Maxwellian velocity distribution of 
the wimps with an rms wimp velocity oi v = 270 km/s. However, there is the possibility of 
a significant uncertainty in the rms wimp velocity with estimates of the uncertainty ranging 
from ±24 km/s to ±70 km/s0. For illustrative purposes we study the effects for the range 
V = 270±50 km/s. The results are shown in Fig.(2a) for the maximum out-going muon flux 
for the Earth and in Fig. (2b) for the maximum out-going muon flux for the Sun. Here we 
find that the ±50 km/s variation of the rms wimp velocity can generate significant variations 



in the out-going muon fluxes. Next we discuss the effect on our analysis of including the 
constraints that will allow the annual modulation signal in the direct detection of wimps 
claimed by the DAMA Collaborationlii. The constraint that the supergravity unifled models 
generate XrV'^oton cross-sections compatible with the annual modulation signal observed 
by DAMA has been analysed elsewhere00 and the reduced set of mSUGRA conflgurations 
compatible with the claimed DAMA signal have been identifled0'0. An analysis including 
the DAMA constraint shows that the peak values of the maximum out-going muon fluxes 
for the Earth and the Sun are not signiflcantly affected. Thus at least in the peak region 
the MACRO limits for the out-going muon flux for the Earth put a further constraint on 
the reduced set of configurations allowed by the DAMA cutS'El. 

In Fig. (3a) we study the mean out-going muon fiux gotten by averaging over the allowed 
parameter space for fixed fine tuning $ as a function of v^- Roughly speaking the fine 
tuning parameter measures how heavy the SUSY spectrum is and we use here the criterion 
given in Ref.0 which defines +](|t- Fig. (3a) shows that the mean out-going muon 

fiux decreases sharply with -\/$- Remarkably one finds that the Earth fiux falls by many 
more decades than the Sun fiux over the allowed range of -\/$- Thus for low values of v^, 
i.e., a/$ < 2.5 (which imply a relatively light SUSY spectrum) the Earth fiux would be 
more easily accessible to neutrino telescope searches, while an opposite situation holds for 
larger -\/$ values, i.e., -\/$ > 2.5 (which imply a relatively heavy SUSY spectrum). This 
means that for a low SUSY mass spectrum, e.g., with values near the current experimental 
lower limits from accelerators, the prime source for the discovery of dark matter wimps for 
neutrino telescopes is the out-going muon fiux from the Earth. However, as the lower limits 
on the SUSY masses at accelerators increase and approach their naturalness limits, the out- 
going muon fiux from the Sun will become the prime source for the discovery of dark matter 
wimps for neutrino telescopes. Of course, as the SUSY spectrum becomes progressively 
heavier, one would need progressively larger sensitivities in neutrino telescopes to probe the 
parameter space remaining in supergravity models. 

Finally we discuss the relative merits of the indirect vs the direct detection (For the 



details of computation of direct detection rates see, e.g., Ref.El). In Fig. (3b) the maximum 
and the minimum out-going muon flux from the Sun and the Earth is plotted against the 
direct detection rate for Ge. The analysis shows that while the direct detection rate spans a 
wide range, i.e., roughly about 4 decades, this range is much smaller than the range of the 
out-going muon flux from the Earth. However, the range of the outgoing muon flux from the 
Sun is comparable to the range of the direct detection rate. The analysis implies that the 
direct detection will certainly be superior to measurements from neutrino telescopes using 
only the out-going muon flux from the Earth at least for relatively larger fine tunings (i.e., 
$ > 2.5, see Fig. (3a)). However, for relatively larger fine tunings the out-going muon fiux 
measurements from the Sun will compete favorably with the direct detection measurements. 

In summary our conclusion is that neutrino telescopes measuring the out-going muon 
fiux only for the Earth or only for the Sun cannot compete with the direct dark matter 
searches in mapping the full parameter space of mSUGRA within the assumed naturalness 
limits. However, the out-going muon fiux measurements from both the Earth and the 
Sun can be competitive with the direct dark matter searches provided the improvement in 
the sensitivities of the indirect and the direct detection measurements keep pace, i.e., the 
sensitivity of measurements improves by the same number of decades in each experiment. We 
have also investigated the effects of non-universalities of the soft SUSY breaking parameters 
at the GUT scale on the analysis presented here. Our conclusions here are very similar to 
the one's given above. The predictions of the muon fluxes given here and their dependence 
on flne tuning will have implications for other neutrino telescopes regarding their ability 
to search for supersymmetric dark matter, e.g., AMANDaEI at the South Pole which is 
collecting data and possibly other future neutrino telescopes which are in the developmental 
stage. 

This research was supported in part by NSF grant PHY-96020274. 
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Fig. 1. (a): Plot of the maximum out-going muon flux for the Earth in mSUGRA as a func- 
tion of the neutralino mass for three different values of the local wimp density corresponding to 
^ = 0.7, 1, 2.3, and jU > (Our sign convention is as in Ref.[19]). The horizontal lines show limits 
from experiment; (b) Same as (a) except for the Sun. 
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Fig. 2. Same as Fig.(la) except that the analysis is for ^ = 1 and for three different rms wimp 
velocities in the halo, i.e., v=220 km/s, 270 km/s, and 320 km/s. (b) Same as (a) except for the 
Sun. 
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Fig. 3. (a): Plot of the mean out-going muon flux for the Sun and the Earth as a function of 
the fine tuning parameter \/$ for three values of ^. The experimental limits from MACRO, which 
are currently the best limits, are shown by the horizontal lines; (b) Plot of the maximum and the 
minimum of the out-going muon flux for the Sun and for the Earth as a function of the direct 
detection rate for Ge. The verticle line is the upper limit from Ref.[20]. 



8 



REFERENCES 



1. V.A. Balkanov et.al., |astro-ph/9903341 



2. M. Ambrosio et.al., |hep-ph/98l2U2(]. 



3. M.M. Boliev et.al., Nucl. Phys. (Proc. Suppl.) 48, 83(1996). 

4. J. Silk, K. Olive, and M. Srednicki, Phys. Rev. Lett. 55, 257(1985); T. Gaisser, G. 
Steigman and S. Tilav, Phys. Rev. D34, 2206(1986); K. Freese, Phys. Lett. B167, 
295(1986); L. Krauss, M. Srednicki and F. Wilczek, Phys. Rev. D33, 2079(1986); 
K. Greist and D. Seckel, Nucl. Phys. B279, 804(1987); S. Ritz and D. Seckel, Nucl. 
Phys.B304, 877(1988): G.F. Giudice and E. Roulet, Nucl. Phys.B316, 429(1989); G.B. 
Gelmini, P. Gondolo and E. Roulet, Nucl. Phys. B351, 623(1991). 

5. M. Kamionkowski, Phys. Rev. D44, 3021(1991); F. Halzen, M. Kamionkowski and T. 
Seltzer, Phys. Rev. D45, 4439(1992); R. Gandhi, J.L. Lopez, D.V. Nanopoulos, K. 
Yuan and A. Zichichi, Phys. Rev. D49, 3691(1994); E. Diehl, G. Kane, C. Kolda and 
J. Wells, Phys. Rev. D52,4223(1995); A. Bottino et.al., Astrop. Phys. 3,65(1995); V. 
Berezinsky et.al., Astrop. Phys. 5, 333(1996). 

6. G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rep. 267, (1995)195. 

7. A.H. Chamseddine, R. Arnowitt, and P. Nath, Phys. Rev. Lett. 49, 970(1982). 

8. For a review see, P. Nath, R. Arnowitt and A.H. Chamseddine, Applied N=l super- 
gravity, Trieste Lectures, 1983(World Scientific, Singapore); H. P. Nilles, Phys. Rep. 
110, 1(1984); H.E. Haber, G.L.Kane. Phys.Rep. 117, 195(1984); R. Arnowitt and P. 
Nath, Proc. of VII J. A. Swieca Summer School, ed. E. Eboli (World Scientific, 1994). 

9. K. Greist and D. Seckel, Nucl. Phys. B283, 681(1987). 

10. P. Nath and R. Arnowitt, Phys. Lett. B336, 395(1994); F. Borzumati, M. Drees, 
and M.M. Nojiri, Phys. Rev.D51, 341(1995); V. Barger and C. Kao, Phys. Rev. 



9 



D57, 3131(1998); H. Baer, M. Brhlik, D. Castano and X. Tata, Phys. Rev. D58, 
015007(1998). 

11. K. Greist and D. Seckel, Phys. Rev. D43, 3191 (1991); P. Gondolo and G. Gelmini, 
Nucl. Phys. B360, 145 (1991); R. Arnowitt and P. Nath, Phys. Lett. B299, 103(1993); 
Phys. Rev. Lett. 70, 3696(1993); Phys. Rev. D54, 2374(1996); H. Baer and M. Brhhk, 
Phys. Rev. D53, 597(1996); V. Barger and C. Kao, Phys. Rev. D57, 3131(1998). 

12. E. Gates, G. Gyuk and M.S. Turner, Phys. Rev. D53, 4138(1996). 

13. G.R.Knapp et.al., Astron. J. 83, 1585(1978); F.J. Kerr and D. Lynden-Bell, Mon. 
Not. R. Astr. Soc. 221, 1023(1986); J.A.R. Caldwell and J.M. Coulsen, Astron. J. 93, 
1090(1987). 



14. M. Kamionkowski and A. Kinkhabwala, Piep-ph/9710337 



15. R. Bernabei et al, Phys. Lett. B424, 195(1998); INFN/AE-98/34, (1998); P. BeUi et.al., 
hep-ph/9gU5mT| . 



16. N. Fornengo,[trep-ph/99045M. 



17. A. Bottino et.al., Phys. Rev. D59, 095004(1999); Astroparticle Phys. 10, 203(1999). 

18. R. Arnowitt and P. Nath, |hep-ph/9902237| (to appear in Phys. Rev. D); M. Brhlik and 
L. Roskowski, |hep-ph/ 9903468 . 



19. K.L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D58, 096004 (1998). For 
other measures of fine tuning see, R. Barbieri and G.F. Giudice, Nucl. Phys. B306, 
63(1988); G.W. Anderson, D.J. Castano and A. Riotto, Phys. Rev. D55, 2950(1997); 
P.H. Chankowski, J. Elhs, M. Olechowski and S. Pokorski, Nucl. Phys.B544, 39(1999). 

20. P. Nath and R. Arnowitt, Phys. Rev. Lett. 74, 4592(1995); R. Arnowitt and P. Nath, 
Phys. Rev.D54, 2374(1996); Phys. Rev. D56, 2820(1997). 

21. M. Beck, (Heidelberg-Moscow Collaboration), Nucl. Phys. Proc. Suppl. 35, 150(1994). 



10 



22. See, F. Halzen, TASI Lectures, |hep-ph/9810368| . 



11 



